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We apply the standard theory of the elastic body to obtain a set of equations describing the 
behavior of an acoustic Gravitational Wave detector, fully taking into account the 3-dimensional 



properties of the mass, the readout and the signal. We show that the advantages given by a Dual 
detector made by two nested oscillators can also be obtained by monitoring two different acoustic 



^ . modes of the same oscillator, thus easing the detector realization. We apply these concepts and by 

H : 

means of an optimization process we derive the main figures of such a single-mass Dual detector 
designed specifically for the frequency interval 2 — 5 kHz. Finally we calculate the SQL sensitivity 
of this detector. 
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I. INTRODUCTION 



Direct detection of gravitational waves (GW), for physics and astrophysics study, is one 
of the great challenges of contemporary experimental physics. Resonant mass bar detectors 



flfl 



the first historically to come to continuous operation, have been improved by 4 
orders of magnitude in energy sensitivity during the last 40 years, so that they can detect 
energy changes of a 2300 kg bar as little as a few hundred of quanta of vibration at about 

1 kHz with a bandwidth up to 100 Hz as. 

Even though the direct and unambiguous detection of GW remains an important goal 
in today's experimental physics, the scientific interest is moving toward the possibility of 
studying the physical and astrophysical features of the radiation. In fact gravitational waves 
will yield unique information about the coherent bulk motions of matter generating the radi- 
ation, revealing features of their source which could not be investigated by electromagnetic, 
cosmic ray, or neutrino studies [6|. In order for acoustic detectors to take part to this "ob- 
servatory phase", a substantial improvement should be achieved both in terms of absolute 
sensitivity and detection bandwidth. In fact only a detector sensitive in a wide frequency 
range, with a good sky coverage for most of the time, will permit to follow the star system 
evolution during a long time interval, improving the signal reconstruction and the detection 
probabilities. Wideband acoustic detectors were recentlyproposed and their expected sensi- 



tivity ana lyti ca lly coated in two idea, configurations flt a precise anaivtica, elation 
of the limit sensitivity of such detectors is still lacking, due to the complexity of the problem. 

Present acoustic gravitational wave detectors consist of a test mass, which is strained 
under the effect of the gravitational wave. A position meter monitors a relevant dimension 
of the test mass, and thereby deduces the equivalent force signal. The sensitivity evaluation 
are traditionally made by modeling the test mass as a simple mechanical oscillator with 
a properly chosen effective mass. This approximation is very good when a narrow band 
detector sensitive around one of its internal resonances is considered: the detector behavior 
within its bandwidth is essentially determined by this resonant mode. On the contrary 
this approximation is seriously limited when applied to a wideband detector. In this case 
the detector is also sensitive at frequencies where many internal resonant modes of the test 
mass will equally contribute to the output signal, each with different amplitudes and phases 
according to their resonant frequencies. The resulting test mass dimension variation is then 
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obtained as the sum of a number of contributions, and can be enhanced or reduced by the 
phase relation of the internal resonant modes. 

Acoustic detectors work at or between the resonant frequencies of the lowest quadrupolar 
normal modes of the test mass, which are set by the geometry and material properties. GW 
sensitive modes in a practical test mass (up to a few meters in diameter) range from about 
1 to several kHz. As a general rule large test masses give better sensitivities and lower 
frequencies, at the expenses of a greater realization difficulty The current predictions of the 
expected rates of GW signals in the kHz band are in rapid evolution, and a full discussion 
of this topic goes beyond the sake of this paper. We point out that a variety of sources have 
been considered recently to give possibly significant GW emission in the acoustic frequency 
range. As an example we give the design guidelines for a detector sensitive in the 2-5 kHz 
range: it would observe a number of interesting phenomena 3, LHJ, [ll| in a frequency range 



121. 



not easily covered by interferometric detectors due to the laser shot noise contribution 

In Section |H] we develop the general mathematical framework of the complete test mass 
+ readout system, leading to a set of formulas in which the response of the full detector is 
related to an external GW force field and to the readout back action force. We also evaluate 
the Standard Quantum Limit (SQL) sensitivity of the detector in the ideal case and propose 
an optimization strategy accounting for the presence of antiresonance frequencies in the 
transfer functions and the noise properties of the readout. Then in section lTTTl we calculate the 
SQL sensitivity of a single-mass Dual acoustic detector specifically designed to be sensitive 
in a broad frequency range. By an in-depth study of the phase relations between the acoustic 
modes, we obtain the geometrical conditions which optimize the detector sensitivity. As our 
calculation takes into account only plane strain resonant modes of the test mass, we validate 
the results by a Finite Element Method (FEM) simulation of the full 3-dimensional body. 
Conclusions and limits of the current work are summarized in Section UVl 



II. WIDEBAND RESONANT GW DETECTOR GENERAL MODEL 

The detector consists of a test mass and a readout system, which monitors a properly 
chosen physical observable of the test mass, namely an average of the surface displacement 
over a sensitive area. Provided that all external noise sources are reduced down to a negligible 
level, the test mass changes its dimensions due to the combined effects of the external GW 
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force, the back action force from the readout and the thermal noise force. In contrast with 
the lumped model case (which holds for a narrow band resonant detector), the effects of these 
forces need to be evaluated differently when the detector is considered as an elastic body. 
In fact the GW equivalent force is applied on the volume of the body, while the back action 
force is applied on the body surface. A complete description of the detector requires the 
knowledge of the position and extension of the surfaces on the test mass which the readout 
samples. In the most general case the readout sensitivity will be position dependent, and a 
weight function must be defined in order to fully describe the interaction with the test mass. 
In the following, vectors are indicated by bold characters. 



A. Test mass mode expansion 

The equations of the motion of an elastic body of density p, forced by a force density 
F(r,t), can be summarized as jlHI ]: 



P^^-LHr,t)] = F(r,t), (1) 



with the appropriate initial and boundary conditions. Here u(r, t) is the displacement field 
of the elastic body and L [u(r, t)} is defined as: 

L[u] = (X + /i) V(V- u) +/iV 2 u. (2) 

The Lame coefficients A and /z depend on the Poisson ratio a v and on the Young modulus 
Y of the material: 

= Yo v = Y 

(l+a p )(l-2a p ) V 2(1 +a p ) ' 1 ' 

We define the displacement normal modes w n (r) as the solutions of the eigenvalue problem: 



-pu 2 n w n = L[w n ] . (4) 

with the boundary conditions on w n given by the requirement that components of the stress 
normal to the test mass surfaces vanish on the test mass surfaces. The normal modes 
constitute an orthogonal complete system, and can be normalized to satisfy the condition: 

/ pw n (r) • w m (r) dV = M 5 nm (5) 
Jv 

where the volume integral is performed on the test mass volume V and M is its total mass. 



The solution of Eq. (P) may then be written as a superposition of the normal modes: 

u ( r '*) = w »( r ) ?»(*)■ ( 6 ) 

n 

The sum runs over an infinite number of elements, as normal modes constitute an infinite 
en— set within the eontinuu, tneehanies —iea. f_ rk Q. The fune- 
tions q n (t) in Eq. © represent the time development of the (n)-th mode, with initial value 
given by: 

0»(°) = T7 / dVpu(r, 0)-w n (r). (7) 

n 

The evaluation of the explicit form of 5 n (i) is straightforward [Ijj] when the driving force 
can be factorized as: 

F(r, t) = G t (t) G r (r) . (8) 

In fact Eq. JIJ becomes: 



and we can apply Eq. multiply by w m and integrate using Eq. (j3J), obtaining: 
M^^- + Mu 2 m q m (t) = G t {t) J dVG r (r)-w m (r). 

The time development of the (m)-th mode is the same of a forced harmonic oscillator. In 
the frequency domain (here and in the following we indicate the Fourier transform with a 
tilda) we have then: 

U") = h i f (UJ \, [ dV G P (r) • w w (r) . (9) 



v 



M (col ~ 

In the hypothesis of dissipation due to the material structure, the system losses can be 



modeled in the frequency domain by including explicitly a damping term 
then: 



16|. We have 



~ ( , 1 G t (u) 



[ #G r (r)-w m (r), (10) 
Jv 



M « - a; 2 ) + iu^u) 
The function 4> n (oj) represents the phase lag between the mass displacement at a given 
frequency and a monochromatic driving force. Experimentally, <p n (u) is found to be roughly 
constant (this is usually referred to as "structural damping" ) and, for low loss materials, it is 
4> n (uj) 1. It can be shown that the normal mode expansion is possible only if the damping 
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term is homogeneous over the test mass volume, as any inhomogeneity of the structure 
datnpin, ca.es a eonpnn g between different .nodes Q. Wnen the .osses ate fte q nene y 
independent the material quality factor Q — 1/0 can be used equivalently in place of the 
phase lag. 



B. Readout 



To detect the external force we measure the resulting strain on some test mass surface S. 
In the small displacement approximation the observable physical quantity, X, of the system 
may be defined as: 

X{t) = [ dsP(r) • u(r,t) . (11) 
Js 

Here P(r) is a weighting function and the integral is performed on the chosen surface S of 
the test mass. We can decide freely the portions of the test mass surface which are sampled 
and their relative weight in the construction of our output variable X(t): the spatial form 
of the weight function P(r) reflects our measurement strategy and detection scheme. For 
example, in an optical readout, P(r) is proportional to the beam spot power profile jl^ . 
while in a capacitive readout [l9| we have P{r) = 1/So on the surface So of the electrodes 
and P(r) = outside (as border effects are usually negligible). 

By Eqs. (Jf3 ITT|) the observable physical quantity of the system is: 

*(*)= y>n(*) / dsP(r) • w n (r), (12) 

or, in the frequency domain: 

*H = %i (13) 

y, [jyW_ G r (r) • w„(r)] [J s dsP(r) ■ w n (r)] 



(ujI - uj 2 ) + iul(f) n (uj) 

This equation allows to evaluate the motion of the system when the specific form of the 
driving force spatial dependence and of the readout weight function are given. Then we 
evaluate the detector response when driven by the relevant forces: a GW, the readout back 
action force and the thermal noise. 

In the case of a linearly polarized gravitational wave propagating along the z-axis, the 



6 



force density applied on the system may be written in cylindrical coordinates (r, 9) as 

G t (t) = \phit) (14) 
G r (r) = W(r) = r cos ( 2 9 + ip ) i r - r sin ( 2 9 + ^ ) i fl , 

where ft(i) is the metric perturbation associated to the GW, ip — for the GW + polarization 
and ip = 7r/4 for the x polarization. When we substitute these in Eq. (jl3j) . the system 
response to a metric perturbation ft becomes: 

X(u) = h(u) H gw (lo) , (15) 

where the detector output is fully described by the transfer function Hcw(u): 

H GW {u) = i (16) 
v [f v dVW(r) ■ w„(r)] [/ s (faP(r) ■ w w (r)] 

The second relevant force applied to the system is the readout back action. Being exerted 
by the readout, this force is applied proportionally to the weight function P(r): 

F BA (t,r) = F BA (t)P(r). (17) 

The test mass response to this input force can again be evaluated by Eq. (JT3j) . if the back 
action surface force Eq. (|T7j) is considered in place of the volume force field G t (t)G r (r): 

X(u) = T ba {u)F B a{u) 

FbM V [/>P(r) ■ w n (r)] 2 

m 2^ K - u?) + MM*) ' 

The transfer function of the system relative to the back action force is then: 

and the resulting displacement can be considered as the coherent sum of the contributions of 
independent harmonic oscillators, each corresponding to a normal mode of the system and 
moving under the same driving force Fba{^)- We note that the displacements contributed 
by different modes could also be in opposite directions, but it is only the resulting sum that 
is physically observable. 



C. Readout optimization and SQL 



The readout noise properties are defined in terms of the single sided power spectral 
densities S xx (u) and Sff(ui). Here S xx (uj) is the equivalent input displacement noise power 
spectral density contributed by the readout; Sff(u) is the equivalent force noise power 
spectrum density due to the readout that drives the system through the transfer function 
Tba [Eq. ([EH) ]- Under the hypothesis of uncorrelated noise sources, the readout contribution 
to the noise on the observable X is: 



Sxx{u) = S xx {u) + \T BA (u)\ 2 S ff {uj) . 



(20) 



A useful figure for expressing the readout performance is its energy resolution expressed 
as number of energy quanta: 

\/S xx {u)S ff (u) 



e r {u) 



h 



The uncertainty relation for a continuous linear measurement requires that 21]: 



e r (u) > 1 



(21) 



(22) 



and the readout is called "quantum limited" when this limit is achieved. It is also helpful 
to define the readout noise stiffness 



K r {u) = yj S ff (u)/S xx (u)) . 



If we now write the power spectral densities S xx {u) and StAuS) in terms of e r and K r : 



the Eq. (|2Uj) becomes: 



S xx (uj) 
S f f{u) 

Sxxiuj) = heJuj) 



1 



+ \T BA {u)\ 2 K r {uj) 



(23) 

(24) 
(25) 

(26) 



It is straightforward to demonstrate that for each u the readout noise contribution is mini- 
mized, giving the SQL [2l( for the detector, if the following equations are both satisfied: 



(27) 
(28) 



Equivalently the following must hold: 




H | T BA (u)\ 

h 



S ff {u) 



T BA {u) b )\ 



(29) 



These relations show that, as obvious, the best detector performances require the use of a 
quantum limited readout and, less trivial, the noise stiffness must be properly matched to 
the test mass mechanical impedance T ba {<jj). The total displacement noise becomes: 



and can be transformed into an equivalent input GW spectral density through the GW 
transfer function H gw {oj) [Eq. (HSJ)]. The detector sensitivity is then: 



We point out that the SQL here derived is calculated taking into account the volume 
properties of test-mass, signal and readout: this is in contrast to the usual derivation of the 
SQL for a GW detector, in which a 1-dimensional model is considered. 

In principle a wide bandwidth detector, designed to be sensitive in a angular frequency 
range uj m i n -ruj max could reach the SQL at every frequency, provided that a readout satisfying 
Eq. (|29~j) on the full interval uj m i n -^-u max can be designed. The elastic body transfer function 
Tba(w) is a rapidly varying function (as we show below in Fig. HJa): sharp magnitude peaks, 
corresponding to the resonant normal modes of the system, alternate with magnitude dips 
called antiresonances. Unfortunately the noise properties of a practical wide bandwidth 
readout have a smooth frequency dependence over the detector sensitivity range and cannot 
match these structures of the transfer function, preventing the possibility of reaching the 
SQL on a wide frequency range. 

The origin of antiresonances is well explained within the modal expansion model. We 
show in Appendix El that in the back action response T ba {uj) there must be an antiresonance 
following a resonance, without exception and regardless of the complexity of structures. An 
optimized readout should have a high back action force noise spectral density at the anti- 
resonances frequencies, where the detector is quite insensitive. On the contrary it must have 
low back action noise at frequencies where resonances of T ba (uj) occur, in order to avoid the 




(30) 




(31) 
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resonance amplification of the readout force noise. The optimization process will approach 
the best performances by averaging among these conflicting requirements. 

The frequency dependence of the readout noise properties depends on its actual design. 
In order to show the optimization strategy, we work out here the simple case of a readout 
with constant noise properties: 

CrM = e (32) 
K r (u) = no- (33) 

The readout contribution to the equivalent input displacement noise on the variable X 
becomes now: 



Sxx{w) = freo 
and the detector sensitivity is simply: 

he 



— + \T ba (u)\ 2 k 

K 



— + \Tba(u)\ 2 Kq 



(34) 



(35) 



\Hgw(u)\ 2 

We find again that for best performances the use of a quantum limited readout (eo = 
1) is required, and the noise stiffness k must be properly matched to the readout. The 
optimization strategy necessarily depends on the characteristic of the expected GW sources. 
In fact for a given GW signal h(u>), the measurements signal to noise ratio S/N is defined 
as: 

s ./ 2 r-iSMi'^. (36) 



N V 71 Jo S hh (uj) 



D. Thermal noise 



In order to reach the sensitivity shown in Eq. (}35j) technical noises must be reduced 
down to a negligible level. Among all possible technical noise sources, the one due to the 
fluctuations of the body volume (the so-called brownian or thermal noise) often ends up as 
a barrier preventing the achievement of the SQL. 

In the case of a system at equilibrium with a bath at temperature T, with a single 
coordinate p and a force F that represents the interaction between the system and the 
externals, the single sided thermal noise power spectrum of the coordinate is predicted by 
the fluctuation-dissipation theorem to be: 

S pp (u) = Sm[r F (w)]. (37) 

10 



where Tf{uj) is the transfer function: 

7> H = m. (38) 
F(w) 

When the system coordinate X is a linear combination weighted by a function P(r), as 
in Eq. the thermal noise on X can be evaluated by the system response to a force 

F(t) P(r) applied to the test mass surface j^]. The thermal noise power spectral density 
on the output variable is then given by the transfer function T ba {oj) [Eq. jUJ)]: 

("> T) = Sm[T B ,( W )] (39) 

a; 

Ak B T v -^M Ij^P(r) ■ w n (r)| 2 
a; ^ M (a£-w 2 ) 2 
The mechanical thermal noise of the system can be reduced in principle below the readout 
contributed noise by reducing the detector temperature T. If the thermal noise evaluated 
in Eq. (J3U|) must be well below the readout SQL noise given by Eq. (|3*0j). then: 

Zm[T BA (u)} « 2h\T BA (uj)\ . (40) 

Notice that the transfer function T b a enters in both sides of Eq. (|40|) . but with different 
meanings: the left-hand side is essentially the incoherent superposition of the thermal noises 
generated by each normal mode of the system, while the right-hand side is the coherent 
superposition of the normal mode responses. 

The effect of the detector thermal noise on the S/N can be easily evaluated if the sensi- 
tivity Shh [Eq. (ESI)] is generalized as: 



S hh (u,T) 



he 


±+\T ba (uj)\ 2 k 


+ S X x{v,T) 







(41) 



where the detector thermal noise Sxx(&,T) is given by Eq. (|39p. In the case of structural 
damping the thermal noise contribution to Shh{^,T) is proportional to the product T(j>. 



III. A CASE STUDY: THE SINGLE-MASS DUAL DETECTOR 

On the basis of the optimization process described in Section 111 CI we can estimate the 
SQL sensitivity of specific detector configurations, provided that the eigenvalue problem Eq. 
(@J) for the test mass is solved. As an example we consider an hollow cylinder test mass (Fig. 
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a) 




Figure 1: (a) Test mass and readout: the displacements induced by a GW propagating along the 
z axis are measured in 8 regions (in black) spanning the whole cylinder height L. (b) Section of the 
detector: the angle a determines the extent of the measuring surfaces and represents a compromise 
between the requirements of wide sampling area and the possibility of implementing 8 independent 
sampling areas on the test mass. The value a = 0.2 rad is used throughout the paper. 

QJi). If we limit our evaluation to GW propagating along the z axis, the symmetry axis of 
the system, the corresponding force does not depend on z: the system response can be well 
described by plane strain solutions, where the displacements are functions of x and y only 
and the displacement along z vanishes. The analytical plain strain solutions are well known 



24j and we use them to evaluate the detector sensitivity; we make use of FEM simulations 



to confirm that the main features of the detector are maintained when a real 3-dimensional 
test mass is considered. 

We show in the following that, thanks to the effect of its lowest order quadrupolar res- 
onant modes, this single test mass detector offers the same advantages of the two nested 
mass configuration (namely the Dual detector), recently proposed ODD : the large band- 
width, set by the frequency difference between the lowest order GW sensitive modes of each 
test mass, and the back action reduction effect, determined by the antiresonance frequency 
placed between these modes. To fully exploit the DUAL detector properties we consider to 
implement a pair of selective readouts. These readouts are sensitive to the two indepen- 
dent components of the quadrupolar modes and reject other classes of modes, which are not 
excited by GWs. 

The selective readouts i? e and i?®(Fig. EI) measures the inner diameters of the hollow 
cylinder and combine them to obtain X 9 = d,2 — d\ and X® = — d%. In each diameter 
evaluation the test mass inner surface displacement is averaged over two opposite areas 
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Figure 2: Sections of the detector showing the detail of the readout measurement strategy. The 
measuring surfaces are arranged in two independent readout configurations: the readout i?® mea- 
sures the difference X^ = d^ — dx, while R® measures the difference X® = d^ — d^. It will be shown 
that R@ and are sensitive respectively to the + and x components of a GW propagating along 
the z-axis. Both readouts are mainly sensitive to the quadrupolar modes of the test mass. 

spanning the full cylinder length. These detection schemes are implemented by the weight 
function: 



where Q(x) represents the unit step function and 2a is the angle subtended by the each 
measured surfaces (Fig. \3p). The value a = 0.2 rad is used in the following: it represents 
a good compromise between the requirements of wide sampling area and the possibility of 
implementing two readouts rotated by 7r/4 on the same test mass. The normalization Sq is 
the area of each measured surface, evaluated as 5*0 — 2ar int L\ as shown in Figure [TJ r int is 
the inner radius of the hollow cylinder. It is straightforward to demonstrate that, according 
to Eq. (|41jl . the detector sensitivity does not depend on the chosen normalization So- 

A. Detector response 

The plain strain normal modes of the motion of the cylinder in the x — y plane are solution 
of the eigenstate problem [Eq. (jlj)], when the following additional boundary conditions are 





(42) 
(43) 



with: 




(44) 
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applied: 



u z =0 

where (w r , tt 2 ) are components of the displacement u(r) in a cylindrical coordinate system. 



The analytical expressions of the norma 
approximation are functions of the kind 



modes of an hollow cylinder in the plane strain 



24): 



w ®n( r ) = fa,n(r) cos (ft 9~) i r + o a , n (r) sin(ft6»)i e (46) 
w ®n( r ) = ~fa,n{r) sin (ft 0) i r + 5f ajn (r)cos(ft6»)i e , 

where the functions f a ,n,9a,n are linear combinations of Bessel functions of the coordinate 
r, with coefficients determined by the boundary conditions. The integer a represents the 
angular symmetry of the mode, while n identifies the mode order within the angular family; 
modes within a given family are ordered by increasing frequency with increasing n. The 
orthogonal displacement fields w® n , w® n represent the same radial distribution of the de- 
formation, mutually rotated by for this reason they share the same eigenvalue uj a , n - We 
call v a>n = uj a>n /(2ii) the resonant frequency of the mode. The displacement w* n is an 
element of an orthogonal complete system, and is normalized to satisfy the condition (with 
s, t = © or <g>): 

/ dVpw s a>n {r)-wl >m {r) =M5 st 6 ab 8 nm . (47) 
Jv 

Any plane strain displacement u p may be written as linear superposition of these basis 
functions, with time dependent coefficients determined by the force acting on the body: 

Up(r,t)= <n(r)<„(t). (48) 

s, a, n 

The plain strain approximation holds in presence of an external force when this vanishes, 
with its first derivatives, along the z— axis: 

G r (r)-i 2 = ^W =0 . (49) 

a z 

Moreover the same conditions must hold for the weight function of the readout used to 
measure the displacement: 

P(r) • i z = = . (50) 

oz 
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The plane strain solutions are exact for an infinite length cylinder, as the fulfillment of plane 
strain condition requires the application, over the terminal sections, of tension or pressure 
adjusted so as to keep constant the length of all longitudinal filaments (that is the body 
height). In the case of a finite length cylinder the solutions can be considered exact only if 
a surface force opposite to the internal stress is applied to the end surfaces of the cylinder. 
According to the Saint- Venant principle Q], we infer that the plane strain approximation 
is good when the length of the cylinder is large compared to its diameter. In this case stress 
and strain in the interior are practically independent of the presence of a force distribution 
applied on the ends, in all the portions of the cylinder except comparatively small slices near 
its ends. 




Figure 3: Modal shape of the first two plane strain quadrupolar modes. The continuous lines 
show the deformation of internal and external surfaces and of 12 radial lines traced ideally inside 
the cylinder section. The shape of the cross section at null deformation is shown as a dashed line. 
Both modes induce deformations of the inner surface of the cylinder and then contribute to the 
detector output, (a) First quadrupolar mode w® 1 (r): both the inner and outer surfaces change 
their shape. The cylinder thickness r ex t — r{ n t remains approximatively constant (fiexural mode), 
(b) Second quadrupolar mode w® 2 (r): the outer surface remains essentially unchanged and the 
deformation involves only inner portions of the cross section. This deformation is then obtained 
by a change of the cylinder thickness (bulk mode) . 

On the basis of the solutions Eq. (|46|). for instance, we can evaluate the effect of a + 
polarized GW propagating along the z-axis, say W + (r), given by Eq. ()14j) with ip = 0. This 
wave can excite only the family of the w® plane strain quadrupolar (a = 2) modes, as we 
have: 

/ dVW+(v).w® n (v) = 0, (51) 
Jv 
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due to the quadrupolar symmetry of the force field. The resulting displacement would be 
the superposition of the displacement contributed by an infinite number of modes with 
a = 2, but we find that the main contribution is given by the two modes of lowest frequency, 
w® 1 (r) and w® 2 (r). These modes induce an elliptical strain of the inner surface of the 
cylinder, as shown in Fig. El We remind that these modal shapes remain unchanged over 
the whole frequency range (namely also far from their resonance), whatever the frequency 
of the driving force, while their relative strength varies according to Eq. (fTU|) . 

As discussed in Appendix |X1 the relative phase between the two contributions depends on 
the sign of the modal constants, which in this case, [according to Eq. (fT7)|) ] are determined 
by the convolution of the GW spatial force W(r) over the modal shape and by the readout 
weight function P e (r) [Eq. (EE21) ]: 

CJ„ oc [ / dVW(r) ■ wf B (r)] [ / dsPM ■ w® n (r)] . (52) 
Jv Js 

The two modes induce a similar deformation of the inner surface, but given by different 

deformations of the cylinder interior: the w® l is a flexural mode and it does not change the 

cylinder thickness r ext — r int , while w® 2 is a bulk mode which changes the thickness. It is 

then not surprising that the analytical evaluation of Eq. (|52~j) results in two modal constants 

with opposite sign. For this reason, when the detector is driven by a GW at frequency v gw 

we recover the characteristics of a Dual detector namely the contribution of the two 

modes: 

- subtracts at low frequencies {y gw < ^2,1), when they follow in phase the driving force. 

- is enhanced when the modes get out of phase in the interval (1^2,1 ^ u 2,2)1 thanks to 
the 7r phase lag of the mode ^2,1 which is driven above its resonant frequency. This 
frequency range determines approximatively the detector bandwidth. 

- subtracts at higher frequencies {y gw > ^2,2) when both modes get out of phase because 
they are driven above their resonance. 

The resulting GW transfer function is shown in Figure 0^: the sensitivity peaks due 
to the quadrupolar modes can be seen clearly. Figure shows the back action transfer 
function T b a, evaluated on the basis of the readout weight function Pe( r ) with a = 0.2 
rad. An antiresonance shows up within the range v 2 ,i ^ ^2,2, as expected on the basis of 
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Figure 4: Transfer functions for the Rq readout with a = 0.2 rad describing the dynamical 
response of an hollow cylinder single-mass DUAL detector (material: molybdenum, outer radius 
r ext = 0.5 m, heigth L = 3 m, form ratio ri n t/r ex t = 0.3). (a) Analytical (continuous line) and 
FEM evaluated (gray line) GW transfer function for a signal propagating along z-axis. Dashed 
line: 3D FEM evaluated response under a + polarized GW traveling along the y-axis. (b) Back 
action transfer function. Continuous line: analytical plain strain evaluation. Gray line: 3D FEM 
simulation. 

the discussion in Appendix |X] The frequency of the antiresonance depends not only on 
the modal constants of the first two quadrupolar modes, but also by the contribution of all 
modes sensed by the readout. We observed that a convergence better than 1% is obtained if 
all modes with frequency below 30 kHz are considered in the evaluation (about 15 modes). 

In Figure 0] we also show the transfer functions obtained by a 3 dimensional FEM har- 
monic a naly sis. These harmonic analysis are performed in the range -j- 10 kHz using the 
ANSYS |2y] full method of solving dynamic equations, proved to be an accurate method 



of evaluating the elastic body dynamical impedance [27[. The FEM code calculates the 
transfer functions from the displacement averaged over the readout surface according to 
the weight function Eq. (|4"2j) . To evaluate Tba, a unit external force is applied onto the 
readout surfaces according to Eq. (JT7J). The GW volume force field Eq. (fTljl is applied to 
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evaluate Hew- As expected, the resonant modes found in the plane strain approximation 
are essentially confirmed by the 3D FEM evaluation, but other modes appear that cannot 
be obtained in the plane strain framework. Thanks to the readout capability of rejecting 
non GW sensitive modes, the displacement induced by these additional modes is efficiently 
averaged out and the detector dynamical performances do not differ from that evaluated in 
the plain strain approximation. This fact confirms that the body dynamics is modified only 
weakly by the absence of the balancing force applied at the bottom and top surfaces of the 
cylinder. 

The FEM code is not limited to plane strain evaluation and can also calculate the sen- 
sitivity pattern of the detector for GWs arriving from an arbitrary direction (9, (p) with an 

I, we separate the 



arbitrary polarization if). Following the approach introduced in ref. 
radiation into its two linear polarization components, described by if) = and if> = 7r/4 in 
the wavefront plane of the GW. When (9, <p) = (0, 0) these force fields reduce respectively 
to the + and x GW described by Eq. (JTJJ) . 

The frequency response to a GW is scaled down for a wave from any direction with respect 
to the z-propagating shown in Figure in the case of a wave traveling along the 

y-axis. We evaluated the scaling factor as a function of the signal direction (9, <f>) and 
polarization if) at a fixed frequency (3000 Hz). For each readout and R®, the normalized 
output for a circularly polarized GW is obtained by quadratically averaging the response 
to if) = and if) = tt/4 waves. The antenna pattern shown in Figure EH refers to circularly 
polarized signals and shows the maximum response between the two readouts. 

The Dual detector is then omnidirectional, in the sense that no blind directions show 
up in its antenna pattern. On the other hand it is not isotropic, and the lowest response 
for a circularly polarized GW is still 30% of the response to an optimally oriented GW, i.e. 
propagating along z. In general, the two linear polarization components of the GW excite 
two independent linear combinations of w® n (r) and wf n (r), which are separately read by 
i?e and R® respectively. Therefore, it is relevant to point out that the single-mass Dual 
detector can resolve both amplitude polarizations of the impinging GW if its propagation 
direction (9, <f>) is known. This is not the case of an interferometric or resonant bar detector, 
which measures only one linear polarization component of the GW. 

The presence of two independent readout channels which monitor different signal polar- 
izations also enhance the possibilities of setting up vetoing strategies in joint observations 
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Figure 5: Antenna pattern of an hollow cylinder single-mass Dual detector (material molybdenum, 
outer radius r ext = 0.5 m, heigth L = 3 m, form ratio ri nt /r ext = 0.3, a=0.2rad) (a) Antenna 
pattern of the detector equipped with the readouts R® and R®. The antenna pattern is represented 
by the modulus of the vector to the plotted surface. Per each direction, we show the maximum 
of the antenna patterns of the readout channels for circularly polarized signals, (b) Section of the 
antenna pattern across the x — y plane. This is the worst case for the detector, which is maximally 
sensitive for GW traveling along the z-axis; nonetheless no blind directions shows up. The response 
ranges from 30% to 55% of the optimal case. 

with other GW detectors, by exploiting Riemann tensor properties to discriminate noise 
events. In addition, the two independent readouts are affected by independent intrinsic 
noise sources, such as thermal and amplifier noises, so that they can be considered as com- 
ing from different detectors as long as environmental disturbances are negligible. 

B. Detector design and optimization 

In this section we determine the optimal detector figures on the basis of the analysis 
developed in Section |TT1 As usual a low dissipation material is required to reduce the effect 
of the thermal noise. Molybdenum represents an interesting choice, as it shows high cross- 
section for GWs and low mechanical losses at low temperatures: ref. |29'| reports a value 
Tcf) < 5 x 10~ 9 K for acoustic modes down to 50 niK. 

In this section we choose to optimize the detector geometry by requiring a specific fre- 
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Figure 6: Frequency of the normal modes in the plane strain approximation. The left axis shows 
the reduced value of the frequency, while the right axis shows the actual value of the frequency 
in the case of a molybdenum test mass with. v ex i =0.5 m and v s = 5650 m/s. The frequency U2,i 
of the first quadrupolar mode (flexural mode) is strongly dependent on the form factor, while 
the frequency ^2,2 of the second quadrupolar mode (bulk mode) is essentially constant. The 
antiresonance frequency (a.r.) is evaluated accounting for the contribution of the number of modes 
needed to obtain a 5% convergence. We point out that while the mode frequencies are solutions 
of the eigenvalue problem Eq. @ and do not depend on the chosen readout, the antiresonance 
frequency is evaluated specifically for the selective readout shown in Figure^. A number of non 
GW sensitive modes are also shown as dotted lines. 

quency band, namely the interval 2-5 kHz. This requirement fixes the frequency of the 
first two quadrupolar modes as u 2) i — 2 kHz and z/ 2j 2 — 5 kHz. As shown in Figure El 
the frequency u 2) 2 depends essentially on the outer radius r ext and on the material sound 
velocity v s . A frequency 1/2,2 ~ 5 kHz could be obtained with a molybdenum test mass 
with r ext = 0.5 m. The height of the cylindrical test mass is fixed to L = 3 m. The inner 
radius, through the form factor ri nt /r ext , determines the frequency of the first quadrupolar 
mode and the placement of the antiresonance within of the bandwidth. In order to better 
understand the effect of the detector geometry on the sensitivity, we evaluate the sensitivity 
of the detector for different values of the form factor in the range 0.1 < r in /r ext < 0.8. 

As discussed in Section III CI the optimization process requires a balancing of the noise 
figures. For this sake we must evaluate as a function of k,q the minimum detectable value 
of a specific GW signal by the use of Eq. (jBTIj) . As in the 2-5 kHz range burst signals are 
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Figure 7: Minimum detectable (at SNR=1) root sum square amplitude h rss of a sine-Gaussian test 
waveform, impinging the detector at 4 different values of the central frequency fo- The calculation 
was performed for a detector form factor ri nt /r ext = 0.3, equipped with a quantum limited readout 
(eo = 1). The value kq ~ 1.7 x 10 11 N/m, indicated by the dotted vertical line, gives the best h rss 
value averaged over the 4 test frequencies. This optimal value of the noise stiffness is independent 
of the T<p value and corresponds to a readout displacement power noise of S xx ~ 6 x 10 -46 m 2 /Hz 
and back action force power noise of Sff ~ 1.8 x 10~ 23 N 2 /Hz. 

expected from compact stellar object, we consider a sine-Gaussian as test waveform: 

(t-tp) 2 

h(t) = h sin[uo (t — to)] e , (53) 

with central frequency f = Uo/^n and r = 2// . For GW detectors it is customary to quote 
performance in terms of the minimum detectable value (at S/N=l) of the root sum square 
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amplitude h. 



TSS ■ 



h 



rss 




dt\h{t)\ 2 




du\h(u)\ 2 . 



(54) 



In order to estimate the detector performances on its full bandwidth, 4 different central 
frequencies were used for the burst of Eq. (J33J). The calculation was also performed as a 
function of the amount of dissipation (T<p value) considered in the test mass. As shown in 
Figure [3 for the case r int /r ext = 0.3 and eo = 1, the optimal value of the noise stiffness is 
k ~ 1.7 x 10 11 N/m. This number is independent of the Tip value and, according to Eq. 
(I24j) . corresponds to a readout displacement power noise of S xx ~ 6 x 10 -46 m 2 /Hz and back 
action force power noise of Sff ~ 1.8 x 10 -23 N 2 /Hz. 

We see in Figure[7|that, at the optimal /t , the thermal noise induced by a finite dissipation 
value Tip = 10~ 8 K reduces the sensitivity of less than 20%. This dissipation is a good 
compromise between sensitivity and feasibility, and can be assumed as the reference Tip 
value for the case eo = 1. In the general case the optimal Tip value is determined by the 
readout energy sensitivity eo- In fact, according to Eq. (|4~T|) . the readout contributed noise 
grows linearly with eo and we find that a value Tip < 10~ 8 e K is needed to limit the thermal 
noise contribution within the 20% of the readout noise. On the contrary the optimal value 
of the noise stiffness k,q is mainly determined by the transfer function Tba{^) and remains 
unchanged. 

Figure |H1 shows the detector performances at the optimal noise stiffness averaged over 
the 4 test frequencies and then plotted against the form factor. Different values of energy 
sensitivity were considered with their corresponding optimal value of Tip. The noise stiffness 
was optimized for each point on the graph and found essentially constant Kq ~ 1.7 x 10 11 
N/m. The dependence on the inner radius is not strong, but the best sensitivity is obtained 
for 0.2 < Ti n /r ext < 0.5, that is when the antiresonance is located approximatively in the 
middle of the bandwidth. The G W sensitivity Shh of this optimal detector (molybdenum, 
Tint = 0.15 m, r ext = 0.5 m, height 3 m) is shown in Figure El 

IV. CONCLUSIONS 

In this paper we discuss wide-band (ie not-resonant) acoustic GW detectors in the 
framework of the elastic theory of continuous bodies and we take fully into account the 
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Figure 8: The detector performances for different values of the readout energy sensitivity is 
averaged over the 4 test frequencies and shown as a function of the form factor ri n /r ext . The error 
bars indicate the scattering of the sensitivity around the average value, a measure of the uniformity 
of the performances over the full bandwidth. A dissipation T(p = 10~ 8 eo K was assumed. 
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Figure 9: Predicted spectral strain sensitivities of the molybdenum detector: weight 22 ton, 
outer diameter 1 m, inner diameter 0.3 m and height 3 m. The curves at different values of 
the readout energy sensitivity eo are evaluated with optimal values of noise stiffness kq = 1.7 x 
10 11 N/m and with dissipation Tcp = 10 _8 eo K. The corresponding displacement noise power of 
the optimal readout is S xx = 6 x 10 -46 eo m 2 /Hz and back action force power noise of Sjf ~ 
1.8 x 1CT 23 e N 2 /Hz. The curve with e = 1 represents the detector SQL. 

3-dimensional properties of signal, test-mass and readout. Once the readout (ie its spatial 
weight function) is chosen, transfer functions are calculated for the GW tidal force, the back 
action noise pressure exerted by the readout and the thermal noise of the test-mass. The 
SQL sensitivity limit of acoustic GW detectors is then derived. 

On the basis of this analysis we design a new kind of acoustic GW wide band detector: the 
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single-mass Dual detector, formed by an hollow cylinder equipped with 2 sets of selective 
readout displaced ir/A apart. This detector would offer the same advantages of a Dual 
detector, namely the kHz-wide frequency range of sensitivity and the back-action reduction 
effect: moreover we show that such a detector has no blind directions for impinging GWs, 
being in this sense omnidirectional even if not isotropic. These findings are analytically 
derived in the plain strain approximation and then extended to the full 3-dimensional body 
by means of FEM numerical analysis. We also show that the single-mass Dual detector can 
measure simultaneously both the amplitude and the polarization of the impinging GW. From 
the practical view-point the single-mass Dual detector removes the technological challenge 
of nesting two independent resonators which is one of the key issues of the recently proposed 
Dual detectors. 

The general treatment of a GW acoustic detector is then applied to the single-mass 
Dual detector in an optimization process that gives the mass dimensions and readout noise 
properties resulting in the best sensitivity in the 2-5 kHz frequency range. We show that a 
molybdenum hollow cylinder single-mass Dual (outer radius r ext = 0.5 m, height L = 3 m, 
form ratio ri nt /r ext = 0.3 and readout angular aperture a = 0.2 rad) can show a sensitivity 
at the level of 1 4- 2 • 10~ 23 / \J Hz in the frequency interval 2 4-5 kHz and truly free from 
noise resonances, provided that T<p < 10~ 8 K, as already achieved for this material. 

The non-resonant readouts for a Dual detector would be evolutions both conceptually and 
in technology of the resonant readouts used in the GW bar detectors: the optomechanical 
one based on Fabry- Perot cavities, and the capacitive one 0], based on SQUID ampli- 
fiers. The main improvements needed to make these readouts suitable in a Dual detector 
configuration can be summarized as: 

The wideband sensitivity. The deformation of the test mass which forms the detector 
needs to be measured by devices that do not show mechanical resonances in the kHz 
wide frequency interval of high sensitivity. 

The selectivity. The Dual detector properties can be fully exploited by a readout scheme 
which is geometrically selective to the fundamental quadrupolar modes. This selec- 
tivity is also effective in reducing down to a negligible level the acoustic modes not 
sensitive to GW that would appear in the detector bandwidth. Moreover the readout 
system must sense the deformation of the resonant masses on a wide surface, in order 
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to be less sensitive to the acoustic modes resonating at higher frequencies, which do 
not carry any gravitational signal. In this way the thermal noise of the detector is 
minimized, while preserving the sensitivity to the signal. 

Without the resonant amplification stage, the displacement sensitivity of readout em- 
ployed in running GW bar detectors ranges in the 10 -19 m/y/Hz decade over a bandwidth 
of hundred of Hz centered at about 1kHz. In order to meet the requirement of a Dual de- 
tector, an improvement is needed by about a factor of 10 4 and these performance must be 
extended to cover the wider Dual bandwidth. For passive readouts, as the capacitor-Squid 
system used in AURIGA, a factor of 10 3 improvement in the displacement sensitivity is 
foreseeable after specific R&D devoted to increase the bias field (presently limited to about 
1% of the material breakdown electric field j^Qj) and to increase the capacitor area (the 
transducer efficiency scales as the square root of the capacitance) and after a reduction of 
the operating temperature down to the 0.1 K range (the SQUID amplifier additive noise is 
linearly proportional to the thermodynamic temperature). For active readouts as the op- 
tomechanical one the 10~ 23 m/ V Hz range is reachable in principle by increasing both 
the finesse (up to 10 6 ) and the light power (up to a few Watts); however such figures could 
end up not to be compatible with a cryogenic operation. Therefore for both technologies we 
can consider as reasonable goal a sensitivity at the level of 10~ 22 m/ \/Hz over the bandwidth 
of a Dual detector. A factor of 10 further could be gained via broadband mechanical lever 
amplifiers, as currently under investigation. 

The wide-area and selectivity requirements could be easily satisfied by capacitive read- 
outs. In fact they are currently implemented on surfaces of 0.05 m 2 , and many capacitors 
could be properly connected to obtain the required total surface and selectivity, as shown in 

n 

ref . j8j . For the optomechanical readout the implementation of selectivity based on crossed 
cavities is under study; in order to meet the requirement of the wide area some of us has 

and its experimental investigation has already started |32j . 
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Figure 10: Typical transfer function for a 4 modes system where the driving force is applied 
with the same spatial weight of the readout. The resonance-antiresonance behavior is due to the 
superposition of the single mode responses and, in the hypotheses of low modal density and small 
damping, an antiresonance is invariably present between the two resonances, (a) Continuous line: 
transfer function obtained when the mode I and mode 2 contributions are only considered. Dashed 
line: transfer function obtained from the contribution of all the 4 modes; the position of the first 
antiresonance is changed by the presence of the higher frequency modes contributions, (b) Transfer 
functions of each mode. 
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Appendix A: MODAL ANALYSIS 

The complex behavior observed in a typical transfer function (Fig. ITUk ) is not surpris- 
ing, as it is well known that structures with low modal density |33| exhibit a pronounced 
resonance-antiresonance behavior in their transfer function 3^ . Antiresonances provide im- 
portant information on the dynamic behavior of the structure under consideration and have 
been applied to several engineering problems, as reported in 35] and Refs. therein. 

Considerable insight may be gained by considering the origin of antiresonances. Within 

26 



the modal expansion model, the system transfer function can be written in the form: 

T{u) = T-2 2 fl n 7 \ > ( A1 ) 

where the explicit form of the modal constants C n and the dissipative factor D n can be 
obtained by a straightforward comparison with Eq. (|13|). As usual the resonant frequen- 
cies represent the zeros of the denominator polynomial, that in the case of small damping 
are approximatively given by the eigenvalues u n . The eigenvalues are determined by the 
mechanical properties of the system and by the boundary conditions applied when solving 
Eq. (J3J). The antiresonance frequencies are the zeros of the numerator polynomial resulting 
from the sum, and their position on the frequency axis depends on the values of the modal 
constants C n . In Figure we show | T(u) | as sum of the contributions of each resonant 
mode for a 4 mode system with Cj > 0. 

Some of the important features of the transfer function can be understood by using a 
simpler example with just two modes and no dissipation: 



TM = ■ (A2) 



If both modal constants are positive, a zero of the numerator polynomial u ar is invariably 
present in the range [cui,^] between the two resonances, at a frequency depending on the 
relative weight between the modal constants C\ and C 2 : 



° 2 (A3) 



Positive modal constants appear in all transfer functions where the driving force is applied 
with the same spatial weight of the readout, as for the back action transfer function Tba [Eq. 
()19j) ]. In this case there must be an antiresonance following a resonance, without exception 
and regardless of the complexity of structures js4|. On the other hand, modal constants 
with opposite signs result in a smooth minimum between the two resonances, located at: 



i 



— + — rr- (A4) 



Modal constants of the detector transfer function to GW [Eq. fllfij) ] can be negative or 
positive, depending on the integrals performed over the modes. In this case the sign is 
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essentially determined by the phase relation between a specific mode and the GW excitation 
force, and by the phase relation between this mode and the weight function. Then we expect 
a mixture of antiresonances and smooth minima. 

We point out that this simple model is only indicative: the actual position of antireso- 
nances and smooth minima between resonances depends on the contribution of the other 
modes, as shown in Figure ITU1 
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